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Abstract 

In this paper, the classical problem of the probabilistic characterization of a random variable is 
re-examined. A random variable is usually described by the probability density function (PDF) or 
by its Fourier transform, namely the characteristic function (CF). The CF can be further expressed 
by a Taylor series involving the moments of the random variable. However, in some circumstances, 
the moments do not exist and the Taylor expansion of the CF is useless. This happens for example 
in the case of a-stable random variables. Here, the problem of representing the CF or the PDF of 
random variables (r.vs) is examined by introducing fractional calculus. Two very remarkable results 
are obtained. Firstly, it is shown that the fractional derivatives of the CF in zero coincide with 
fractional moments. This is true also in case of CF not derivable in zero (like the CF of a-stable 
r.vs). Moreover, it is shown that the CF may be represented by a generalized Taylor expansion 
involving fractional moments. The generalized Taylor series proposed is also able to represent the 
PDF in a perfect dual representation to that in terms of CF. The PDF representation in terms 
of fractional moments is especially accurate in the tails and this is very important in engineering 
problems, like estimating structural safety. 

1 Introduction 

In many cases of engineering interest, it is useful to describe strength or mechanical properties of materials, 
geometrical features of structures or loads and so on, as random variables. A random variable is fully charac- 
terized by the PDF or by its spectral counterpart, namely the CF, or by the moments (or cumulants) of every 
order. Moments and cumulants are related to the coefHcients of the CityplaceTaylor expansion in zero of the 
CF and of the log CF, respectively. 

Yet, the moment representation of the CF is not always feasible. Indeed, if the derivatives in zero of the 
CF do not exist, the Taylor moments series is meaningless. For example, the CF of a-stable random variable 
[1], [2], is not derivable and consequently such representation does not exist. Another example is the response 
to non-linear structures under parametric stochastic input, which may have divergent moments starting from 
a certain order, even if the system is stable in probability. With such moment structure, the CF cannot be 
restored. From these observations, it might be concluded that moment resolution strategies are often infeasible. 
In this paper, we investigate on the particular class of moments with complex exponent and we show how 
powerful this extended class is. To the author's knowledge, such moments have not been investigated in 
literature. Then, definitions and existence conditions of such moments will be given in the following. Of 
course, as the PDF is a real function, it follows straightforwardly that if the moments up to a given integer 
order n exists then all the complex moments with real part less than n and greater than zero must exist. The 
question is: are such moments useful in reconstructing the CF and the PDF? 

In order to answer to this question fractional calculus will be adopted. The latter has received a growing interest 
in the last century and applications in physics and biophysics, in quantum mechanics, in the study of porous 
systems (gathered in the book [3]), in fracture mechanics [4], in non-local elasticity [5], [6], to cite just few, are 
available in literature. In stochastic dynamics, the PDF of response to differential equations driven by Levy 
a-stable white noise processes is ruled by a fractional differential equation, involving fractional derivative in the 



*Publication info: Cottono G., Di Paola M., On the use of fractional calculus for the probabilistic characterization 
of random variables, Probabilistic Engineering Mechanics, Volume 24, Issue 3, 2009, Pages 321-330, ISSN 0266-8920, 
10.1016/j.probengmech.2008.08.002. 

tE-mail: giulio.cottone@tum.de; giuliocottone@yahoo.it 



1 



diffusive term [7]. Fractional derivatives are encountered also in random vibration with frequency dependent 
parameter [8] and in the analysis of linear or non-linear systems driven by fractional Brownian motion [9], [10]. 
By means of the fractional calculus, we show that the fractional derivatives of the CF in zero coincide with 
fractional moments, also in the case of CF not derivable in zero (like the CF of a-stable r.vs). To the authors' 
knowledge, the only author giving a relation between fractional moments and fractional derivatives of the CF 
in zero is Wolfe [11] by using fractional derivatives. However, the expression obtained by Wolfe is given in an 
integral form quite different from the classical expression relating moments and derivatives in zero of the CF. 
Indicating with C the set of complex numbers, we will show that fractional derivatives and integrals of order 
7 € C of the CF calculated in zero are nothing more than particular moments of order 7 of the random 
variable. Once this remarkable result is achieved, the usefulness of fractional moments for the probabilistic 
characterization of random variables is pointed out by using a generalized Taylor integral theorem proposed by 
Samko et al. [12], which involves fractional derivatives. It is shown that, by means of a suitable discretization, 
the series here proposed involving moments of complex order looks like the classical Taylor series. A satisfactory 
representation of the CF is made possible also for the cases in which the CF exhibits a slope discontinuity in 
zero as in the case of a-stable random variables. Moreover, it is well known that, when the classical Taylor 
series is truncated, the CF exhibits unsatisfactory trends on tails, whilst in this paper it is shown that the 
fractional Taylor series always fulfils the desired fundamental property that the CF vanishes at infinity. By 
some easy algebra we will show a very useful dual representation of the PDF by means of complex moments. 
Finally, we will show by numerical examples that with a finite set of complex moments one can have very good 
approximations both of the PDF and of the CF. In particular, the PDF is very good approximated in the tails 
and this is important in problems like the estimation of structural safety. 

2 Preliminary Concepts and Definitions 

In this section, some well-known concepts on the probabilistic characterization of random variables, as well as 
some definition of fractional differential calculus, are briefly summarized for clarity's sake and with the aim to 
introduce appropriate symbols. 

Let G R be a real random variable whose probabilistic characterization may be given both by the PDF 
Px ix) and by its Fourier transform, namely the CF (f)x {'&), that is 



(j)x = E [exp [i-dX)] = / exp [Mx) px {x) dx (1) 
./ — 00 

where i? £ R, i = -\/— T is the imaginary unit and E [•] indicates average. Provided that moments E [X-*] with 
j — 1, 2, defined as 

E \X'] = I px ix)x'dx (2) 
exist, then (j)x (i^) can be expanded in Taylor series 

'^-w = E^^M (3) 



J'- 



due to the property 

d'(tix W 



E 



From equation (U), it is clear that E [X-*] exists if the j-th derivative in zero of the CF exists. For example, 
the moments E ^X\^\ of the a-stable random variables do not exist for p > a and, since the stability index a 
ranges from zero up to 2, then the Taylor expansion ((3} cannot be applied. In spite of this, since E [\X\^] and 
E [X^] £ C with p < a exist, is it possible expanding the CF into some series involving fractional moments of 
order p < a? The answer is affirmative, as we prove after recalling few remarks on fractional calculus. For more 
exhaustive treatment on fractional calculus, readers are referred to the excellent encyclopaedic book of [12], 
and to [13], [14]. The greatest difficulty one has to overcome dealing with fractional calculus is represented by 
so many definitions of fractional derivatives present in literature, and each definition has its own peculiarities 
and technicalities. As it will be stated clearly in the next sections, every definition is good for our purpose. In 
particular, we recall the definitions of the Riemann-Liouville, Marchaud and Riesz fractional derivatives and 
integrals since they are useful for the ensuing results. 

The Riemann-Liouville (RL) fractional integrals of order p € M. > 0, denoted as {l±f) (x) are defined as 
follows 

(/i/) (x) -i-^ r e-'f T dC (5) 



r(p) Jo 

where F (p) /|^°° ^''^^exp (— ^) d^ is the Euler gamma function that interpolates the factorial function, that 
is (n — 1) \ — r(7i). Often, the operators (x) and {Itf) (x) are referred as left and right hand sided. 
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respectively. In some literature, these integrals {l±f) (x) are also denoted as Liouville-Weyl (LW) fractional 
integrals. The RL fractional derivatives, denoted as (x), are expressed as 

with p £ R > and n = [p] + 1, being [p] the integer part of p. Examining RL definitions it is clear that the 
conditions of existence of such operators depend strongly on the behavior of the function / (x) at ±oo. 
The Marchaud definition of the fractional derivative, denoted as (D^-/) (x) is expressed as 

fD" f) (x) - r 

^^■^^^^-rii-{p})l, e+M 

where /'Ip'^ (a;) denotes the derivative of order equal to the integer part of the real number p, and {p} = p—[p]. 
It is worth to note that the Marchaud fractional derivative exists also for functions growing at infinity as 
\x\P-\e>0. 

The Riesz fractional integration, denoted as (/''/) (x) is defined as follows 

2r(p)cos(pV2)i_|e-x|^-'' 
with p G K > 0, p 7^ 1, 3, 5. ...The Riesz integral may be expressed by the RL operator as follows 

and the Riesz fractional derivative denoted as (C/) (x) may be represented in terms of Marchaud fractional 
derivative [12] 

iV^f) (-) = -J^^^) ((DU) (-) + (D^f) (x)) (10) 

Moreover, the Riesz fractional derivative may be also expressed in terms of RL fractional derivatives as 
follows [9], [10] 

i^^f) (-) = - 2eos(pV2) + ^''^ 

The operation of fractional integration, which has been introduced for real p > 0, can be made meaningful also 
for complex value of the index, say it 7 = p + 177, with p,?) G R. Indeed, the fractional integrals {l±f) (x) and 
{P f) {x) are properly defined under the condition that p > 0. In the same way, fractional derivatives (V^-f) (x), 
{V^ f) (x) are defined accordingly to relations ([71, (fTIJ)l . (fTTj) . having care that n = [Re^f] = [p]. It must be 
clear that integrals or derivatives of complex order 7 G C (and p 0) represent an analytic continuation in the 
parameter 7 of fractional integrals and derivatives originally defined for rj — Imy — 0. Sufficient conditions on 
the existence of such operators and extensions to the case of Re"/ — are reported in Samko et al.([12], p. 38 - 
39). Here, we stress only that the fractional derivatives are derivatives of convolution integral that produces a 
smoothing on the original function: then, the fractional derivative of a function might exist even if the function 
is not classically differentiable. 

A very important feature of fractional operators is their behaviour with respect to the Fourier transform. 
The Fourier transform of a function / (x) defined in the whole real axis, is a function of the variable i9 defined as 
T{f (x) ■ ■&} = e'^'^f (x) dx, jointly with the inverse transform J^-i {g {j}) ■ x] = (27r)-^ e"^'"'^ g (i9) d-&. 
Then, it has be proven that the Fourier transform of Riemann-Liouville, Marchaud, and Riesz fractional 
derivatives of order 7 are given, respectively, by the relations 

T { {VI f) (x) ; ^} = (Ti^y T {/ (x) ; ,?} (12) 

T { (Dl/) {x) ■ ^} = (Ti^)" T {/ (x) ; ^} (13) 

T{{V-^j) (x) ; ^} = - 1^1^ (x);^} (14) 

where Re'y > (see [12], pp 137; 218). 

Analogously, the Fourier transforms of Riemann-Liouville and Riesz fractional integrals are given by the 
relations 

T {Ilf {x) ■ ^} = (Ti^)-' T{f{x);^},0< Rej < 1 (15) 
T{rf (x) ■ ^} = \^rT{f (x) ; ^} , < Re^ < 1 (16) 

Some of the main properties and rules of the aforementioned fractional derivatives and integrals like Mellin 
transform and composition rules are reported in Appendix A. 
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3 Fractional derivatives and integrals of the characteristic func- 
tion 

In this section, useful relationships between fractional moments and Riesz fractional integrals and derivatives 
are derived. In the following of the paper, it will be indicated by 7 a complex number, ^ — p + i-q, with p, 77 
reals. The function used in the following {^^ix^ has to be understood ([12], p. 137) as 

{Tix)"' = exp ^7 In Ixj =F -^'Sgn (x)^ , 7 £ C, Re-y > 0, x (17) 

being 

C 1 X > 

sgn (x) =■'1 X = (18) 
[ -1 2; < 

the usual signum function of a real variable. Now, following ea. H12l) . the Fourier transform of the Riesz fractional 
derivative of the CF is given as 

T { {Vl<f>x) {^) ; x) = {-^ixy F{<l,x (1?) ; x} , i?e7 > (19) 

Applying the inverse Fourier transform to both sides of ea. H19|) gives 

/oc 
{Tixy px{x)eyip{i'dx)dx (20) 
- OC 

Finally, by evaluating eq. (|20[) in -i? = we obtain the fundamental relation 

/oo 
{^ixy px (x) dx = {Vl<t)x) (0) ; 7 e C, Rej > (21) 
-oo 

stating that the fractional RL derivative of order 7 of the CF in zero equals the complex moment of order 7 
of the random variable X. That is, complex moments are ruled by an expression very similar to that obtained 
by classical moments (see eq.(|4]). The easy procedure leading to (|21l) may also be used for demonstrating the 
following identities, starting from the Fourier properties p2l) - (ll6p 

(O>A-)(0) = -£[|XP],i?e7>0 (22) 

ir<t)x) (0) = E [jXr^] , i?e7 > (23) 

{-Dl<Px) (0) = E [iTiXy] , i?e7 > (24) 

(40x) (0) = E [{TiXyy ,Re^>0 (25) 

We want to remark that, Wolfe [11] derived fractional moments, by using the Marchaud fractional derivatives 
in the form 

= [(D'+<^x) (^)|,=, - {-Dlcbx) (^)U=J ^} + (26) 

+ i-''Re{{Dl^x) WU^o} 

with — cxj < q < 00. Notice that the structure of the ea. H26|) hides the simplicity of the conceptual connection 
between moments and characteristic function present in the corresponding ea. (|22|l . Furthermore, ea. (|26|l does 
not give fractional moments explicitly, because of the integral involved. Eas. (|21|l - (|25|l have been obtained in 
the most general context by using complex exponents and can be written for real exponent as well. 

It should be noted that moments of the type E [(=piX)*^] , after some cumbersome algebra, may be ex- 
pressed in terms of moments E [X^y , in the form 

E [{iXy] = 2 cos ) £ [\X\y - ^-^E [X^] (27) 



E [i-iXy] = i-''E [X''] (28) 

E [iixyy = 2cos E [\xyy - pe [x-y (29) 
E [(-ixyy = i'' E [x-y (30) 

with Re-y > 0. 

Summing up, in this first step we have shown that the classical relation between CF and moments (|4]) is 
extended also to complex moments, by the RL definition of fractional operators. The next step is to use a 
representation of the CF (and of the PDF) in terms of its fractional derivatives of complex order calculated in 
zero. It will be demonstrated that fractional moments can be used in reconstructing the CF (or the PDF) of 
every distribution (also distributions with discontinuity in zero in the CF) extending the well known eq.(|3}. 



4 



4 Generalized Taylor series using fractional moments 



There are many generalizations of the Taylor series involving fractional derivatives, see e.g. [12], [15]-[21]. 
Riemann himself made the first attempt. In his posthumous published work he wrote 

r m + r 

/(- + '^)= E r(^ + . + i) (31) 

m— — oo 

with r a fixed real number and {'D2^ f) (x) is the left-handed RL fractional derivative on a finite support 

If r in eg. pifl is an integer, all the terms m < — r disappear and it is equivalent to the classical Taylor 
series. However, Hardy [17] proved that the expansion H31|) is a divergent series. For h large enough the 
series converges, but, as we are interested on the values of the derivatives of the CF in zero, the generalized 
Taylor expansion of Riemann is useless for our scope. Many other generalizations of the Taylor expansion using 
fractional derivatives exist. In Appendix B some of them are analyzed in detail showing that they have some 
undesired property for the representation of the CF. 

In this section an integral form of Taylor expansion related to the inverse Mellin transform, proposed by 
Samko et al.([12], pp. 144-145) will be applied to the CF and to the PDF of a random variable. Then, by means 
of eq. (|2ip - (|25|l . generalization of eq.Q will be provided. 

The starting point is to recognize that the RL integral (O may be interpreted as the Mellin transforms (see 
Appendix A) of the functions / (a; ± ^), with x fixed, that is 

p oo 

r (7) (^)= / r"\f(^±c)d^ (33) 

Under this perspective, by performing the inverse Mellin transform, one gets 

1 /-p+ioo 

f{^±i)-7^. r(7)(J^/)(x)|C|-^d7,P = i?e7>0 (34) 

^ TT I J p-ioo 

that for a; = becomes 

1 /'p+ioo 

f{±0 = 7^. r(7)(/;/)(o)|c|-^d7 (35) 

As claimed in [12], this integral may be interpreted as an integral form of the Taylor expansion in the sense 
that knowing {llf) (0) on some line p = Re [7] > 0, then the value of the function / in the whole axis may be 
predicted. Under this consideration, ea. (|35p is the generalized Taylor expansion proposed by Samko et al. 

5 Integral representation of the CF in terms of complex mo- 
ments 

Let us suppose that (^Il^cfix) (0), with 7 G C and p — Re"f > 0, exists. In the following, conditions of existence 
will be given. Then, particularizing ea. (|35|l and using ea. (|25|) for the CF one obtains 
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p+i 00 



'Px (T^) = 7r- r (7) (/J0x) (0) i^r^ dj = (36) 

1 rp+ioc 

= r{^)E[{±rX)-'']\^r'dr, 

■^"^ J p-ioo 

The line integral in eq. (|36[l can be rewritten as 

(T^) = ^ T{p + i^)E [(iiX)-^-'"] l^r''-"' dn (37) 

Some remarks are necessary at this point: (i) due to the position p > 0, the moments E |^(=FiX)^''^*''j in 
ea. p7[l have intrinsically negative real order; (ii) random variables with divergent integer moments do have, at 
least in an open interval, negative real order moments and can be represented by ea. ()37|) : (iii) ea. (|36|l can be 
evaluated for every value of p inside an interval, that is < p < 1, as we prove in the following. 

For clarity's sake, we develop an easy example by considering the CF of a Gaussian distribution showing the 
relation between the generalized integral representation (|36|l proposed and the Taylor series Q. We recall that 
the characteristic function of a standard Gaussian random variable X is (fix {'&) ~ e""* ^'^ . In order to apply 
ea. p6p the value of p = Re"f must be properly selected. As first consideration, from the definition of the RL 
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fractional integral, the operator (/J^x) {"&) is meaningful with Re"/ > 0. In addition, in case of Gaussian r.v., 
moments of the type E [{'^iX) '] — J^^px (x) {^^ix) ' dx exist in the range — oo < Re"/ < 1. Combining 
the two conditions of existence, one obtains the so called fundamental strip < Re^f < 1, that is an open 
interval in which ea. (|36p holds. It will be shown that every distribution can be represented by the integral 
representation p6l if Rej belongs to the range < Re^f < 1. Further, the integrand function is holomorph 
inside the fundamental strip, ensuring that one can choose every value inside the fundamental strip. Outside 
the fundamental strip the integrand is not necessarily holomorph and shows isolated singularities in the no- 
positive part of the real axes. In Figure 1, just the real part of the integrand is plotted choosing a particular 
value of 1? = I? > 0. 

It has to be noted that, in this particular case, the function has poles at 0,-2,-4,-6, ...etc. Now, due to the 
residue theorem, integrating along the line with real abscissa < p < 1, corresponds to integrate in the whole 
half plane at the left of the line, that is to the sum of the residues at the poles of the integrand, obtaining: 

<bxW = — r (7) E [i-^X)-^] d-'d-i = (38) 

— oo 

Res {V (k) E [(-jX)-'=j ■d-''^ = 

k-0,-2,... 




where Res(-) is the residue, and the property Res(r(— fc)) — {—l)^/k\, k = 0,1,2, ...has been used. The 
former expression is the connection between the Taylor series representation of the characteristic function and 
the integral representation proposed. 

Now, it is well-known that using the Taylor series in terms of integer moments one needs a closure scheme 
in order to properly truncate the sum to a finite number of terms. In the integral representation by means of 
complex moments it is easier to find in which interval the integrand does not contribute anymore to the value 
of the integral. Indeed, as Figure 1 highlights, the integral p7l can be approximated by truncating it in a range 
[—fj, ??] and then numerically evaluated. In the next section, performing a basic rectangle integration scheme it 
is shown that with a finite number of fractional moments one can approximate either the CF both the PDF of 
every distribution. 



5.1 Series approximation of the CF in terms of fractional moments 

The integral we want to evaluate numerically by means of a simple integration scheme is 

<^x(±i?) = ^ I^^T{p + zr^)E[{TiX)-''-^'^]m-'-^'^drj (39) 

that, can be approximated by truncating the extremes 

0A- (±1?) = 2V /''_r(p + i77)iJ [(T^X)-"-"'] \&\-^-'^dr^ (40) 

Then, setting rj = kA, with fc G N and A G and letting 7^ = p + ikA, with p > 0, and being fj = mA (see 
Figure 2, ea. (|40p may be rewritten in discrete approximated form as follows 

A 

<?ix (1?) = — r ilk) E [i^iXy-"'] W--"- fon9 > (41) 

fc— — m 

A 

(1?) ?^ — Y ^ ^ l^r^' fori? < (42) 

Of course, CF has the property that {—'&) is the complex conjugate of (i?), in the following denoted 
as <l)x {1}), and this fact will lead to the interesting simplification that only one equation between eas. (|4ip - (|42|l 
suffices to represent the whole CF. 

By introducing eqs. (|29|l - (|30p in (|4ip - (|42|l . the series representation of the CF in terms of moments and 
absolute moments may be rewritten in the form 

0x (i?) ~ V r (7fc) f'-E [X-'^] l^]-'" (43) 

fc— — m 
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for i9 > 



A 

-^^W-^ E r(7.){2cos(^)i5[|^r"''] + (44) 

fc — — OO 

-i'O'E [X-''"] l-ff]-^"} 

fon? < 

Note that, since 4}x ~ (t>*x i^'^), ea. (|43|l is enough in order to define the CF in all the range ]— oo, oo[. 
Furthermore, ea. (|43|l has a mathematical form very similar to that of the classical Taylor expansion of the CF 
given in eg.©. 

We investigate now for the bounds of p such that this series representation is valid. Of course, due to the 
presence of the RL fractional integral, p must be positive from definition ([S]), so, the first lower bound of p is 
zero. In order to have a higher bound, one shall recall the condition of existence of the inverse Mellin transform. 
The choice of the value of p, indeed, must belong to the so called fundamental strip of the Mellin transform in 
ea. p6|l . (see Appendix A). In order to define the fundamental strip of the Mellin transform, it suffices to note 
that the integral 

poo 

M {^x (±0 C'^x (±0 (45) 

converges, if the integral 

< r e-'\<t>x{±0\d^ (46) 
Jo 

at the r.h.s. converges. As the characteristic function is absolutely convergent in the real axis, ea. (|46p is 
bounded at least for every < p < 1. Then, in order to use the integrals (|39|) - H40|) or the series ()41|) - H42|) . it 
suffices to choose a value of p in the interval < p < 1. This represents the strictest condition of existence 
of the series proposed. Of course, there are many distributions satisfying H45p in a wider fundamental strip. 
For example, the CF of a standard Gaussian distribution admits a Mellin transform in the fundamental strip 
< p < oo and we might choose every p > in evaluating eQ. (|39|) . 

Suppose we want to use eg. 1)43 [I for the representation of the CF of an Q-stable random variable. It is 
known that for such distributions moments of the type E [X"''^"'] are finite, although complexes, only if 
— 1 < — p < a, otherwise they diverge due to the heavy tails of the PDF. Being the stability index a defined in 
the range < a < 2, the condition < p < 1 on the fundamental strip allows to represent the CF by means of 
such moments E [X~''~^^] with intrinsically negative real order —1 < — p < 0. 

In order to understand the role played by the parameter A, we first recall that the series form representation 
is the numerical integration of ea. (|36p . For fixed p, the integration is performed in the imaginary axis as shown 
in Figure 2. Small values of A produces high accuracy, however since from a numerical point of view in any 
cases the summation (|43|1 has to be truncated retaining a certain number, say m number of terms, if A is very 
small then mA = fj remains small and we may exclude in the integration significant part of the integrand. 
Then, fj must be chosen in a such way that F (p + if]) E (=piX)~''~''' is negligible. 



Jo 



5.2 Probability density function representation in terms of fractional mo- 
ments 

In engineering problems, i.e. in solution of stochastic differential equations or in path integral methods, we 
are concerned more in probability, rather than in its spectral counterpart. In this section it is derived a dual 
integral representation of the PDF in terms of complex moments, by means of Fourier transform. 

The inverse Fourier transform of the CF can be performed by {4>x {"&) ',x} or, that is the same, 
J-~^ {(jix (i? > 0) + 4>x (i? < 0) ; x} by the linearity of the integration involved. Then, the inverse Fourier trans- 
form of relation H36|l may be performed in the following way 

(47) 

1 rp^ioo f POO 

px (x) = {0x W ; x} = — / r (7) (/;0x) (0) / exp i-i^x) 

(2-K) I J p-ioo L Jo 

+ {Il(t>x) (0) j {-^y exp i-i-dx) di?! d-y 
that, under the condition < p < 1 can be further simplified 

P^(^) = 7T4r-/ r(7)r(l-7){(/>x)(0)(ia;)--^ + (/Z<^x)(0)Hx)--^} d7 (48) 
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Finally, substituting eq. (|25p in the latter equation, the relation searched reads 

(^) = / r (7) r (1 - 7) {E [{-^X)-''] [ixy-' + E {-^xy-'} dj (49) 

and it is the integral form of the Taylor approximation of the PDF in terms of moments. 

It has to be remarked that, since (f)x = (l>*x then px (x) may be also rewritten in a more compact 

form as follows 

PX {x) = —Re / r (7) r (1 - 7) (/>x) (0) {ixy-^ d^ = (50) 

f p-\-i 00 



1 r rp+ioa 

= -±j-Re / r (7) r (1 - 7) £ [{-ixy^] [ixy-' d^ 

Then, with the same set of moments E [(— iX)^^] one can represent both the CF, with ea. ()36|l evaluated only 
for J? > 0, or the PDF evaluated by ea.(|50[l. without additional calculation (i.e. Fourier transform). 
Expression (|5U|I can be approximated in the form 

Px{x)^^Rel f2 r{j,)r{l~j,)E[{^iX)-''']{ixy^A (51) 

K. k— — 77i J 

Taking into account ea. (|3Q|) . previous expression simplifies 

Pxix)?^^Rel r(7fc)r(l-7fe)i^'=i5[X-^^] (jx)^'^-4 (52) 

V fc— — m J 

Such a perfect duality between the two representations in i9 and x domains, evident comparing eqs. H43|) and 
(|52|l . may not be evicted working in terms of classical moments. In fact an inverse Fourier transform of (|3]) 
gives the PDF px (x) in the following form 

that cannot be if practical use in representing the PDF because of the Dirac's delta 5 (x) derivatives. By 
summing up knowledge of moments of the type E [X "=] , 7^ G C are able to represent both the CF and the 
PDF in the whole 1? and x domains. 



6 Numerical results 

In this section, we present some simple applications in order to show the validity of ea. (|2H) for both the cases 
in which the CF is differentiable in zero and for case of a-stable random variables. In the following, p will be 
used to indicate the real part of 7. As it has been already pointed out, the value of p inside the fundamental 
strip can be arbitrarily chosen, then the results here reported for particular values of p can be reproduced of 
course with different choice of its value. Let X be a random variable with uniform distribution in [—a, a], then 
all the moments exist and are given as 

2k 

E [X^*^] = £ [X''+'] = 0, fc = 1, 2, ... (54) 

while the CF is 

'da 



sin (i?a) . 



The classical Taylor expansion of the CF is then 



fc=0 



Such a series converges uniformly to (fix (1?) in the entire domain. However, when a truncation is performed 
the CF diverges as it is shown in Figure 3a. The latter circumstance produces an undesirable result when 
the PDF of X has to be restored by Fourier transform of the truncated series H56I) . The converse happens 
by using the fractional series proposed. In fact, as shown in Figure 3b, by applying (|41|) - H42|) to the CF 
(|55|l . one notes that the fractional series remains always convergent in the whole domain and for i9 — >■ 00, 
(fx ("!?) — > 0. Easy calculations give the moments of complex order (—7) for the uniform distribution, i.e. 
E [(±iX)^^] = (ia) ' /(I — 7) defined in a strip —00 < p < 1; in particular, in Figure 3b it has been assumed 
that p = 0.4 and A = 0.4. 
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In order to show that the proposed series works also for complex CF, in Figure 4, the proposed expansion 
of the CF is reported for a Rayleigh distribution with scale [da vedere]. In Figure 4a and 4b the real part 
and the imaginary part are plotted, respectively. The complex moments of a Rayleigh distribution with scale 
parameter ct G R > are given by E [(ipiX)"^] = 2~^^^ (rpia)"^ T (1 — 7/2), within the strip — cxj < p < 2. 
Then, in this case one could choose the fundamental strip in the interval < p < 2. Despite, the knowledge of 
fractional moments in some line inside the strip < p < 1 ensures that the series (|4ip - (l42|l is meaningful. The 
results of the fractional series reported in Figure 4 have been calculated with the choice of p = 0.4, but similar 
results may be obtained also for other values in the range of validity of < p < 1, and A = 0.4. 

A more interesting and striking example is the symmetric Cauchy random variable X with px {^) = 
TT^^ /{x^ + 1). For such a distribution the moments exist only for — 1 < p < 1, then the usual Taylor expansion 
fails because it would be just a sum of divergent terms, and therefore meaningless. Further, in this case it may 
be shown that moments E [(=piX)~^] = 1, in the strip —1 < p < 1, and the approach in this paper proposed 
using eqs. (|4ip and (|42|l . is a very good approximation, as shown in Figure 5. 

Then, we remark that the fractional series representation works also for distributions with moments existing 
only in a limited range, like the stable distributions, because one needs only the existence of a real moment 
E [X^''] in, at least, one real value p, such that < p < 1. 

To enforce this concept, as last example, the Levy stable variable with PDF given by 
Px exp {—l/2x), whose moments E [{—iX) ''] exist if the condition —1/2 < p < cxj 

is satisfied. The CF of this distribution has been approximated with p = 0.9, then moments of the type 
E [X""'^-*''] , have been used in series (|43l)-(|44} or, that is the same, moments E [(=FiX)"°"""'] in (l4T|) -- (|42ll . 



Figure 6a and Figure 6b show the convergence of the series proposed. 

Now, we report some interesting applications of the fractional series ea. (|51|) to find the PDF from the 
knowledge of some moments, that is, without using inverse Fourier transform. In Figure 7 one can see that 
the PDF of a Gaussian random variable is well approximated by 30 complex moments in the whole domain, 
having used p — 0.4 and A = 0.4. Figure 8 and Figure 9 show that in order to represent the PDF of a Cauchy 
and a Levy random variables respectively, 10 complex moments are sufficient, with the choice of p = 0.4 and 
A = 0.4. 

7 Conclusions 

In this paper, a new perspective of the probabilistic characterization of the random variables has been high- 
lighted. It has been shown that the fractional derivatives of the characteristic function in zero coincide with the 
fractional moments. Then, integral Taylor approximation due to Samko et al. [12] has been the starting point 
to develop a representation of both the density and the characteristic function of every random variable. Such 
a new representation ofi'ers many advantages: (i) it may be applied to not-differentiable in zero characteristic 
functions, and to densities with heavy-tails; (ii) with a unique set of moments E [{—iX) '] both the PDF 
and the CF can be represented in dual form; (iii) the series proposed are not of integer power-type and once 
truncated they vanish at infinity. 

The authors have recently extended to multivariate random variables the representation proposed in this 
paper. Moreover, a path integral method for the solution of non-linear one-dimensional stochastic differential 
equation excited by white noise processes based on complex moments has been studied (see [22] and [23]), 
while the solution of multidimensional stochastic differential equation is underway by the authors. Indeed, the 
path integral method, which is based on the evolution in time of the response density, has proven to take full 
advantage from the PDF representation in terms of complex moments presented in this paper. 

8 Appendix A 

We gather in this appendix some well-known definitions of the Mellin transform needed in the text. Further, 
we briefly describe the so called composition rules of fractional operators. 

The Mellin transform ([12], p. 25; [24], p. 41) of a function f{x),x > 0, s G C, is defined as = 
M{f{x);s} = !S° fiO C'rf? along the inverse operator /(a;) = M W {s) ■ x} = {2niy^ ip (s) x-'ds, 

with c — Re (s). The Mellin transform if (s) of the function / (x), such that 



exists, and is analytic, in the fundamental strip —p < Re [s] < —q. In order to perform the inverse Mellin 
transform c = -Res must be chosen inside the fundamental strip. Tables of Mellin transforms of commonly used 
functions are given in [24] and [25]. 

Properties of nested application of the fractional operators, called composition rules, are useful in calcula- 
tions. Suppose the existences of the integrals and derivatives involved in the following relations, and consider 




(57) 
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Re'yi > and Re'y2 > 0. Then, the properties 

J71 ^72 J ^ J71 +72 J ^ J72 J71 J (gg) 
Jll Jl2 J _ j7lH-72 J _ J72J71 y (59) 

on the composition of different order integrals, are valid. On the converse, composition rule between 
fractional derivative and fractional integral does not follow straightforwardly. Indeed, to better understand 
the commutation property between fractional integration and differentiation some parallel with the standard 
differential calculus might help. In ordinary calculus, it is well-known that performing firstly the integral of a 
function and then a derivative produces the original function, i.e. ^ / (x) dx = f {x). The same happens 
to n-order derivatives and n-folded integrals. Exactly in the same way, applying the fractional operators in 
the order 'VI''/ gives the original function /. 

Instead, performing firstly the derivative of / and then integrating, i.e. ■£:f{x)dx, one obtains the 
function plus a constant, or for n-order derivatives and n-folded integrals one obtains the function / plus a 
polynomial of order n — I. This simple argument is valid also in case of fractional operators and consequently 
the relation /''C/ — f is not in general true, but depends on the function / (x). It has been shown ([12], p. 
43-45) that the latter is true only for those functions vanishing at the interval boundaries. 



9 Appendix B 



Although many generalization of the Taylor series exist (see e.g. [15] - [21]), the only series adapt to express 
the CF in terms of moments is the series proposed, descending from the Mellin transform. In order to avoid 
useless attempts to readers interested on this topic, we describe in this short appendix some reasoning on other 
possible ways, at least at first sight, to develop a CF representation by means of a generalized Taylor expansion. 

The first family of generalized Taylor series we examine descend from the composition rule of fractional 
operators aforementioned. In order to simplify the discussion we will use the left-handed RL fractional integral 
and derivative on a finite support defined as 



It has been already stressed that, in general, the fractional derivatives and integrals do not commute, if the 
derivative is performed as first operation, that is 

l2+V:^f7^fix),Re^>0 (62) 

but, it has been proved in [12], that (|62|) is replaced by 

j _ -,7-fc-l 

I2^v:j = / (x) - ^ r(7-fc) D-'-'ll-^''^ (a) , (63) 

where the operator D-' is the classical derivative of order jth. Rearranging this relation one find an analogue of 
the Taylor theorem, expressing the function / (x) as sum of fractional derivatives calculated in a point 

f = E r(7 + fc + i) " + I'^r^^Tf, (64) 

with n — [Rej] + 1. From this formula one can state that, if for a function the property I2_^1^2^f = / i^) holds, 
then the function cannot be represented by the fractional Taylor series (|64p . because all the terms of the sum 
actually vanish. In the perspective of this paper, in order to understand if relation (|64|l is useful, it suffices to 
study the operation I2^T>2_^4>x {'&)■ Recalling that the CF of a random variable X is (jix (t?) = E [exp (Mx)], or, 
that is the same, IJ^D^+^x (i^) = II+T^a^E [exp {i'dx)]=E [/2+^a+ exp (i'dx)] . By means of direct calculations 
it is easy to verify that Iq^Dq^ exp{Mx) = exp{i-dx) and therefore series (|64|l is useless to our scope, with 
a = 0. Other values of a would not produce a connection between CF and moments. Same reasoning applies 
to the generalized Taylor series proposed by Dzherbashyan and Nersesyan [19]. 
Other Taylor expansions have been given by Jumarie [20] 



/ (^) - E r(pf+i) (0) ' < P < 1 (65) 

or, for z £ C, by Osier [16] 

fi^)-pj: Fp^(</)M. (66) 
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Using relations (|65p . (|66|) it was not possible to extend the eq. Q for two reasons. Firstly, it was not possible to 
find an explicit relation between the fractional derivatives ^Dg+l^-f ) (^) of the CF and the fractional moments, 
being ^Dq^(^x^ defined just for values of > 0, since the Fourier transform of such operator, like p2p . does 

not exist to our knowledge. Moreover, even if such an expression would exist, terms of the kind (t>q'^(I>x^ {'&), 
would involve fractional moments whose real part increases with k = — 1, 0, 1, ... and therefore not applicable 
to CF not difi'erentiable in zero. 
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Figure 1: Plot of the real part of the integrand in ea. (|5S)) for ^ = ^ > 0. 
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Figure 2: Numerical integration scheme in the fundamental strip. 
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Figure 5: Fractional Taylor series of the CF of a Cauchy distribution with zero location and unitary 
scale parameter. 




Figure 6: Fractional Taylor series of the CF of a Levy distribution with zero location and unitary scale 
parameter. 
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Figure 7: Fractional Taylor series of the PDF of a Gaussian distribution with mean 2 and unitary 
standard deviation. 




Figure 8: Fractional Taylor series of the PDF of a Cauchy distribution with zero location and unitary 
scale parameter. 
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Figure 9: Fractional Taylor series of the PDF of a Levy distribution with zero location and unitary 
scale parameter. 
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